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We study physics potential of Very Long Base-Line (VLBL) Neutrino-Oscillation Experiments with the High 
Intensity Proton Accelerator, which will be completed by the year 2007 in Tokai- village, Japan. As a target, a 
100 kton-level water- Cerenkov detector is considered at 2,100 km away. Assuming the pulsed narrow-band 
beams, we study sensitivity of such experiments to the neutrino mass hierarchy, the mass-squared differences, one 
CP phase and three angles of the lepton-flavor-mixing matrix. We find that experiments at a distance 2,100 km 
can determine the neutrino mass hierarchy if the mixing matrix element |f7 e 3| is not too small. The CP phase 
and | C/ e 3 1 can be constrained if the large-mixing-angle solution of the solar-neutrino deficit is realized. 



1. Introduction 

In order to measure the neutrino oscillation 
parameters, such as the neutrino mass-squared 
differences and the elements of the 3x3 Maki- 
Nakagawa-Sakata (MNS) lepton flavor-mixing 
matrix elements jlj , various long base-line (LBL) 
neutrino oscillation experiments are proposed. In 
Japan, as a sequel to the K2K experiment, a new 
LBL neutrino oscillation experiment between the 
High Intensity Proton Accelerator (HIPA) |2| and 
the Super-Kamiokande (SK) with the base-line 
length of L=295 km has been proposed ||. The 
HIPA has a 50 GeV proton accelerator, which will 
be completed by the year 2007 in the site of Japan 
Atomic Energy Research Institute (JAERI) as 
the joint project of KEK and JAERI. It is ex- 
pected to deliver 10 21 protons on target (POT) 
in one year. The intensity of the neutrino beam 
in the ~ 1 GeV range is two orders of magni- 
tude higher than that of the KEK beam for the 
K2K experiment. The HIPA-to-SK experiment 
with (E u )= 1.3 GeV will measure the larger mass- 
squared difference with 3 % accuracy and the 
mixing angle at about 1 % accuracy. 

In this report we examine physics potential of 
Very Long Base-Line (VLBL) neutrino oscillation 
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experiments with HIPA and a huge neutrino de- 
tector j|] in Beijing, at about L=2,100 km away. 
As a beam option, the pulsed narrow-band v 
beam (NBB) is assumed. For a target at L=2,100 
km, we consider a 100 kton water- Cerenkov de- 
tector which is capable of measuring the z/„ to 
v e transition probability, Pv^—tu^i an d the v „ sur- 
vival probability, P v _, v . We study the sensi- 
tivity of such experiments to the neutrino mass 
hierarchy, the mass-squared differences, the three 
angles and one CP phase of the MNS matrix || . 

2. Neutrino oscillation in the three- 
neutrino model 

The MNS matrix has three mixing angles and 
three phases in general. We adopt the following 
parameterization 
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Two Majorana phases in the latter matrix, ip 2 
and ip 3 , do not contribute to the neutrino oscilla- 
tion. The former matrix can be parameterized by 
three mixing angles and one phase just the same 
as the CKM matrix. Because the present neutrino 
oscillation experiments constrain directly the el- 
ements, U e 2, t/ e 3, and these three matrix 
elements are taken as the independent param- 
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eters. Without losing generality, we can take 
U e 2 and to be real and non-negative. By 
allowing U e 3 to have the complex phase U & 3 — 
|C^e3|e _,s,5 MNs (0 < S MNS < 2ir) , the four indepen- 
dent parameters are U e 2, U^, \ U e 3\ and <5 MNS . All 
the other matrix elements are then determined by 
the unitary conditions. 

The constraints on the MNS matrix elements 
and the mass-squared differences are given by 
the atmospheric- and solar-neutrino and the 
CHOOZ reactor experiments. An analysis of the 
atmospheric-neutrino data from the SK exper- 
iment §] finds sin 2 26» ATM - (0.88 - 1.0) and 
Sm 2 ATM (eV 2 ) ~ (1.6 - 4.0) x 10" 3 . From the ob- 
servations of the solar-neutrino deficit by the SK 
collaboration @, the MSW larg e-mixing-angle so- 
lution (LMA) is preferred to MSW small-mixing- 
angle solution (SMA), MSW low-c$m 2 solution 
(LOW) and Vacuum Oscillation solution (VO). 
The CHOOZ experiment || gives the constraint 
sin 2 26 CHOOZ < 0.1 for <5 TO 2 HOQz > 1.0 x 10" 3 eV 2 . 
From the above experiments, the independent pa- 
rameters in the MNS matrix are obtained as 
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Since all the above constraints are obtained 
from the survival probabilities which are even- 
functions of 5m 2 j, there are four neutrino-mass 
hierarchy cases corresponding to the sign of the 
5m 2 j as shown in Fig. [[} If the MSW effect is rele- 
vant for the solar-neutrino oscillation, then the hi- 
erarchy cases II and IV are not favored. The hier- 
archy I (III) is called 'normal' ('inverted') hierar- 
chy, which corresponds to 8m\ 2 > and 8m\ 3 > 
( (5m 2 3 < 0). From the Super-Nova 1987A ob- 
servation, the hierarchy I is favored against the 
hierarchy III 0. 



Figure 1. Schematical view of the four cases of 
neutrino-mass hierarchy. 



3. VLBL neutrino oscillation experiment 
with L=2,100 km 

We explore the capability of a VLBL experi- 
ment with L=2,100 km by using u M NBB. The 
parameterization of the NBB is made to study 
the effects of changing the peak energy, E p , of the 
NBB. We find that the use of two different-energy 
NBB's, NBB's with E p = 4 GeV (NBB(4GeV)) 
and E p = 6 GeV (NBB(6GeV)), improves the 
physics resolving power of the experiment signif- 
icantly. Fig. H shows our parameterizations of 
the main part of NBB(4GeV) and NBB(6GeV) 
fluxes (x neutrino-energy). Also overlayed in 
Fig. H are the oscillation probabilities, which 
are calculated for sin 2 28 ATM — 1.0, <5"i 2 TM = 

2 „;„2 nQ X™2 10 X 

0.1, <5 MNS = 270° and 
the matter density, p — 3 g/cm 3 , with the hierar- 
chy I. Note that the NBB(6GeV) makes the e-like 
event rates high while keeping /i-like event rates 
low. 

For a detector at L=2,100 km, we assume a 
100 kton water-Cerenkov detector. It has a capa- 
bility of distinguishing e charged-current (CC) 
events from CC events, but does not distin- 
guish their charges. The detection efficiencies of 
the /i-like and e-like events are assumed to be 100 
% for simplicity. The statistical errors of each 
predictions are then simply square values of the 
observed numbers of events. We do not require 
capability of the detector to reconstruct the neu- 
trino energy. 



3.5 x 10~ 3 eV 2 , sin 2 20 SOL = 0.8, <5m 2 OL 
10" 5 eV 2 , sin 2 20 CHOO5 
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Figure 2. Overlayed are the NBB's (flux times 
E v ) used in our analysis and the neutrino oscilla- 
tion probabilities at L=2,100 km. 



The signals of our analysis are then the num- 
bers of Up CC events and those of v e CC events 
from the NBB. For NBB(^ p ), they are calcu- 
lated as 

N(j i {e),E p ) = 

MN A [ " dE v <f>(E p )aff e) P v ^„ Me) , 
Jo 

where M is the statistical significance (mass of 
the detector (= 100 kton) times the running 
years), Na = 6.0 1 7 x 10 23 is the Avogadro num- 
ber, is the flux at L = 2, 100 km. The 
* s ^ e u v>( Ve ) ^ cross section per nucleon 
off water target. 

Fig. U shows the numbers of signal events, 
N(fi, E p ) and N(e, E p ), expected at the base-line 
length of L=2,100 km for 500 kton-year with (a) 
the NBB(4GeV) and with (b) the NBB(6GeV). 
The predictions for the four cases of the neutrino- 
mass hierarchies, I, II, III and IV, are indicated 
explicitly. We take 

sin 2 20 ATM = 1.0, Sm 2 ATM = 3.5 x 10~ 3 eV 2 
sin 2 20 CHOOZ = 0.06 ,0.1, 
6 MNS - 0° - 360° , 

with p = 3 g/cm 3 . The remaining two pa- 
rameters, sin 2 26* SOL and <5«t 2 OL , are taken as 
the following sets for the three possible solu- 
tions to the solar-neutrino deficit anomaly : 



(sin 2 20 SOL ,<5m 2 OL ) = (0.8,15(3) x 10~ 5 eV 2 ) for 
LMA, (7 x 10" 3 , 5 x 10" 6 eV 2 ) for SMA, (0.7, 7 x 
10~ n eV 2 ) for VO. 

For the hierarchy I, the predictions reside in 
the corner at small N(fj,,E p ) and large N(e,E p ). 
For the LMA scenario, the predictions are shown 
by circles when <5 MNS is allowed to vary freely. 
We show the four representative phase values 
by solid-circle (S Mtqs — 0°), solid-square (90°), 
open-circle (180°), and open-square (270°) in the 
largest circle. The two grand circles with smaller 
N(/j,,E p ) give the predictions of the LMA solu- 
tion for Sm^ OL — 15 x 10~ 5 eV 2 , and those with 
larger N(fj,,E p ) are for <5m 2 OL = 3 x 10~ 5 eV 2 . 
The v e CC event number N(e,E p ) grows with 
increasing sin 2 26* CHOOZ . It is clearly seen that 
<5 MNS dependence is larger for larger <5w 2 QL and 
for larger sin 2 2# CHOOZ . Note that the S MNS de- 
pendence of N(e, 4GeV) is 'orthogonal' to that of 
iV(e, 6GeV). The predictions of the SMA and the 
VO parameters appear just above the upper LMA 
circles, where the <5 MNS dependence diminishes to 
zero. We cannot distinguish the VO predictions 
between the neutrino-mass hierarchy I and II, nor 
between III and IV. 

All the predictions of the hierarchy cases III 
and IV have very small N(e,E p ), typically a fac- 
tor of 5 or more smaller in magnitude than the 
corresponding number of events in case I. This 
striking sensitivity of P u on the hierarchy 
cases is the bases of the capability of distinguish- 
ing the cases in the VLBL experiments. On the 
other hand, we find that the 5 % level differences 
in N(/j,, E p ) between the hierarchy cases I and III 
are not useful for this purpose because P Vfi ^ iU ^ 
depends strongly on the atmospheric parameters 
5m 2 TM and sin 2 26> ATM , and also because N(/j,, E p ) 
suffers from the neutrino-beam flux uncertainty. 

We take into account for the background con- 
tributions to the /i-like and the e-like events from 
the secondary-beams and the r pure-leptonic de- 
cays. The e-like events also receive contributions 
from the neutral-current (NC) events where pro- 
duced 7r°'s mimic electron shower in the water- 
Cerenkov detector. The error in the e/NC 
misidentification probability, P e /NC = 0.0055 ± 
0.00055, is accounted for as a systematic error. 
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(b) NBB(6GeV) 
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Figure 3. The expected signal numbers, N(fi,E p ) and N(e,E p ), at L=2,100 km for 500 kton-year with 
(a) NBB(4GeV) and (b) NBB(6GeV). The predictions for the four-cases of the neutrino-mass hierarchies 
are depicted as I, II, III and IV. The hierarchy cases II and IV are shown for completeness sake even 
though they do not lead to the MSW solution of the solar- neutrino deficit. 



The v r CC events with hadronic r-decays are 
counted as NC-like events with the same -P e /NC 
factor. The 10 % errors in the branching fractions 
of the r decay are accounted for as systematic er- 
rors. 

In addition, we account for the uncertainties in 
the total fluxes of the neutrino NBB's and the 
matter density as the major part of the system- 
atic uncertainty. We allocate common 3 % uncer- 
tainty for all the NBB's and 3.3 % overall uncer- 
tainty in the matter density. 

4. Results of the \ 2 analysis 

We examine the capability of the VLBL exper- 
iments with L — 2, 100 km in determining the 
model parameters. The following questions are 
of our concern. 

1. Can we distinguish the neutrino-mass hierar- 
chy cases I (Smf 3 > 0) and III (<5m 2 3 < 0) ? 

2. Can we measure sin 2 26> CHOOZ and r5 MNS ? 

3. Can we distinguish the solar-neutrino oscilla- 
tion scenarios ? 

4. How much can we improve the measurements 



We examine the questions by combining two 
experiments with different base-line length, L = 
2,100 km (HIPA-to-Beijing) and L = 295 km 
(HIPA-to-SK). For a VLBL experiment at L = 
2, 100 km, we assume the statistical significance 
of 500 kton-year each with the NBB(4GeV)and 
NBB(6GeV). As for the LBL experiment at L= 
295 km, we assume that 100 kton-year for the low- 
energy NBB with (p w >= 2 GeV (NBB(2tt)). The 
data obtained from the LBL experiment is what 
SK can gather in approximately 5 years with 10 21 
POT par year. (The latest design of the HIPA- 
to-SK project Q adopts off-axis beams in the 
first stage.) 

The x 2 is a function of the three angles, the two 
mass-squared differences, the CP phase, the flux 
normalization factors and the matter density. For 
a given set of these parameters, we have defined 
predictions of the /i-like and e-like event numbers 
including both the signal and the background. 

4.1. I v.s. Ill 

We show in Fig. ^ the minimum % 2 as a func- 
tion of sin 2 2# CHOOZ by assuming the hierarchy III 
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in the analysis, when the expected event numbers 
are generated for the following input (true) pa- 
rameters with the hierarchy I in the LMA region. 

Input (true) parameters : 
sin 2 26 true = 1.0, 5m 2 true = 3.5 x 1CT 3 eV 2 

ATM ' ATM 

sin 2 2(9*™ e = 0.8, Sm 2 J T rue = 10 x 10~ 5 eV 2 
sin 2 2^™ e oz = 0.02,0.04,0.1, 
§ true =o°, 90°, 180°, 270°. 

MNS ' 3 ' 

8m\ 2 = (Sm 2 ^™ 6 , Sm 2 3 — <5m 2 T '™ e (hierarchy I) 

In total 12 minimum x 2 lines are labeled by 
the input values of sin 2 20*™* , = 0.02 (thin 

CHOOZ * 

lines), 0.04 (medium-thick lines), 0.1 (thick 
lines) and = 0° (solid lines), 90° (long- 

dashed lines), 180° (short-dashed lines), 270° 
(dot-dashed lines). 



Hierarchy I (true) vs III (fit) 




0.02 0.04 0.06 0.08 0.1 0.12 



sin 29 CHOOZ 



Figure 4. Minimum \ 2 as a function of 
sin 2 2(9 CHOOZ when the hierarchy III is postulated 
in the analysis. The event numbers are calculated 
for the hierarchy I in the LMA region. The input 
(true) values of sin 2 CHOOZ and <5 MNS are shown 
by the line thickness and the line type, respec- 
tively. 



The x 2 fit has been performed by assuming the 
hierarchy III. The minimum of the \ 2 function is 



found for a sin 2 2# CHOOZ value in the range below 
0.12, by varying the parameters, sin 2 20 SOL and 
rjm 2 OL , in the LMA region, and the remaining 
three parameters, sin 2 20 ATM , £m ATM and <5 MNS , 
freely. 

Fitting parameters: 

sin 2 20 ATM : free, Sm 2 ATM : free, 

sin 2 26> SOL = 0.7-0.9, <5m 2 OL = (3-15) x 10~ 5 eV 2 , 

sin 2 26» CHOOZ : free, <5 MNS : free, 

Sm 2 2 = ^ m soL'^ m i3 = — ^ m ATM (hierarchy HI) 

The uncertainties in the total fluxes of the NBB's 
and the matter density are taken into account. 

The minimum X 2 for sin 2 26»^ H " e oz = 0.1 is al- 
ways larger than 50, which allows us to distin- 
guish the hierarchy I from III at 7cr level. We 
also find that the hierarchy III can be rejected 
at 3rj level when sin 2 26* CHOOZ > 0.04, and at lu 
level sin 2 20 CHOOZ > 0.02. The capability of dis- 
tinguishing the hierarchy cases I and III comes 
from the remarkable differences of the transition 
probability P v —>v e between them. 

4.2. sin 2 26» CHOOZ and S MNS 

Fig. shows the minimum \ 2 in the 
sin 2 2# CHOOZ v.s. <5 MNS plane. The event num- 
bers are calculated at sin 2 20* r "^ = 0.06 for 

CHOO Z 

<^™ e = 0°, 90°, 180°, 270° in each figure. The 
other parameters are taken as the same as in 
Fig. ^. We show the input parameter point by 
a solid-circle and the minimum x 2 = 1,4, and 
9 by the solid-line, dashcd-line and dotted-line, 
respectively. The x 2 fit has been performed by 
assuming the LMA scenario with the hierarchy I, 
while the other parameters are freely varied. 

When <5^™ e is around 0° or 180°, 0.04 < 
sin 2 2# CHOOZ < 0.1 is obtained and rJ MNS can 
be constrained to local values at la level. On 
the other hand, when S tru ^ is around 90° or 

5 MNS 

270°, 0.02 < sin 2 26> CHOOZ < 0.08 or 0.045 < 
sin 2 2# CHOOZ < 0.12 are obtained, respectively, 
but we cannot constrain <5 MNS . As shown in Fig. 
the S MNS = 90° (270°) point for the hierarchy I 
lies in the lower (higher) N(e,E p ) corner. The 
same values of A(e,4GeV) and iV(e,6GeV) can 
be obtained for the other 5 MNS by changing ap- 
propriately the sin 2 2(9 CHOOZ parameter. On the 
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5. Summary 

The HIPA will be completed by the year 2007 
and deliver 10 21 POT at 50 GeV for one year 
operation. Assuming the NBB's from the HIPA 
and a 100 kton-level water-Cerenkov detector, we 
study the physics potential of the VLBL experi- 
ments with L=2,100 km as a sequel to the pro- 
posed HIPA-to-SK LBL experiment with L =295 
km. Thanks to the enhancement of matter effect, 
it is expected to distinguish the neutrino-mass hi- 
erarchy cases in such VLBL experiments. We find 
that a combination of LBL experiments at L—295 
km and VLBL experiments at L=2,100 km can 
determine the neutrino mass hierarchy at 3a level 
if sin 2 26> CHOOZ >0.04. Also it is found that, if the 
LMA scenario is realized in Nature, sin 2 2# CHOOZ 
and 5 mNS can be constrained at la level in favor- 
able cases. 



Figure 5. Minimum \ 2 contour plot in the 
sin 2 26» CHOOz v.s. <5 MNS plane. The four values 
of the input parameters are given by the solid- 
circles. The \ 2 fit nas been performed by assum- 
ing the LMA scenario with the hierarchy I. 



other hand, in the case of 5 MNS = 0° or 180°, 
it is difficult to reproduce both N(e, 4GeV) and 
N(e, 6GeV) simultaneously because of the or- 
thogonal S MNS -dependence between V(e,4GeV) 
and 7V(e,6GeV). 

4.3. Solar- and atmospheric-neutrino pa- 
rameters 

We give only results for the remaining two ques- 
tions. If the LMA scenario is realized in Nature 
and sin 2 26^™ e oz > 0.06, the SMA/LOW/VO sce- 
narios can be rejected at la level when S tru ^ is 

J MNS 

around 0° or 180° but not at all when S true is 

MNS 

around 90° or 270°. 

For the atmospheric-neutrino oscillation pa- 
rameters, sin 2 26* ATM is measured to 1% level and 
<5m 2 TM with the 3% accuracy when sin 2 20*™^ = 
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